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ABSTRACT 

Moir6  interferometry  with  line  densities  of  1200  and  40  lines  per  mm  was  used  to 
determine  the  two  orthogonal  displacements  surrounding  a  stably  extending  crack  in  a 
2024-T3  aluminum  alloy,  single  edge  cracked  specimen.  The  test  protocol  consisted  of 
using  the  fine  moird  grating  prior  to  and  up  to  the  onset  of  crack  extension  and  the  coarse 
moir9  grating  for  the  ensuing  crack  extension  up  to  Aa  •  6  mm.  The  displacement  fields 
were  used  to  compute  the  J-integrals  for  various  contours  during  crack  tip  blunting  and 
crack  extension.  As  expected,  the  far-field  J-integral  value  prior  to  stable  crack  growth 
coincided  with  the  LEFM  strain  energy  release  rate,  G,  and  va’  fated  the  experimental 
procedure.  However,  the  J  values  obtained  from  the  near  tip  contour  increased  slowly 
while  the  far  field  J  values  increased  rapidly  with  increasing  stable  crack  growth.  The 
HRR  displacement  field  was  computed  from  the  experimentally  determined  far  field  J. 

The  HRR  displacement  field  agreed  with  the  measured  displacement  field  prior  to  stable 
crack  growth  since  J  -  G.  However,  the  HRR  horizontal  displacement  field  progressively 
deviated  from  the  measured  values  with  crack  extension. 

KEYWORDS 

Elastic-plastic  fracture  mechanics,  J-integral,  HRR  field,  SEN  fracture  specimen, 
moird  interferometry. 

INTRODUCTION 

For  the  past  five  years,  the  authors  and  their  colleagues  have  used  experimentally 
determined  displacement  fields  to  compute  directly  the  J-integral  in  thin  aluminum 
fracture  specimens  (Kang  and  Kobayashi,  1988;  Dadkhah  and  Kobayashi,  1990,  1992,  1993: 
May  et  al,  1993).  The  contour  integration  was  performed  using  the  definition  of  the  J- 
integral  (Rice,  1968)  with  an  added  assumption  that  the  elastic-plastic  response  of  the 
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specimen  material  could  be  represented  by  a  power  hardening  form  with  experimentally 
determined  power  hardening  coefficients.  These  J-integral  values  agreed  well  with  the 
corresponding  elastic  values  under  low  load  and  with  the  known  elastic-plastic  solutions 
(Kumar  et  al,  1981)  at  higher  loads  prior  to  stable  crack  growth.  Similar  conclusions 
were  obtained  by  Schultheisz  (1991)  who  studied  the  effect  of  crack  tunnelling  in  lO-mm 
thick,  4340  alloy  steel,  three  point  bend  specimens.  Under  a  small  stable  crack  growth  of 
1  to  3  mm,  however,  the  near-field  J-integral  values  increasingly  deviated  from  the 
known  solutions  (Kumar  et  al,  1981)  and  were  as  much  as  one  quarter  of  that  computed  by 
the  aforementioned  EPRI  Elastic- Plastic  Handbook  values. 

These  J-integral  values  were  then  used  to  compute  the  HRR  displacement  field 
(Hutchinson,  1968:  Rice  and  Rosengren,  1968)  for  comparison  with  the  experimentally 
measured  displacement.  Theoretically  the  HRR  displacement  field  should  have  coincided 
with  the  measured  displacement  field  within  the  experimental  and  numerical  error 
bounds.  In  reality,  the  measured  and  HRR  displacements  perpendicular  to  the  crack,  i.e. 
the  v-displacements,  agreed  well  but  substantial  differences  were  found  in  the 
corresponding  values  for  the  displacements  parallel  to  the  crack,  i.e.  the  u-displacements. 
Attempts  were  also  made  to  characterize  the  second  order  term  in  the  crack  tip 
asymptotic  displacement  field  which  is  defined  here  as  the  difference  between  the 
measured  and  HRR  displacements.  These  results  suggested  that  neither  the  two  parameter 
characterization  by  the  Q-stress  (O'Dowd  and  Shih;  1991,  1992,  1993),  nor  the  T-stress 
(Betegon  and  Hancock;  1991),  prevailed.  Thus,  without  the  presence  of  an  HRR  field  or  an 
HRR  field  with  the  higher  order  terms  at  the  elastic-plastic  crack  tip,  the  much  heralded 
J-integral  loses  its  physical  significance  as  the  strength  of  the  HRR  singular  field. 

Many  theoretical  and  numerical  studies  related  to  the  J-integral  and  the  HRR  field 
have  been  published  since  1968.  Much  of  these  studies  were  limited  to  computing  the  J- 
integral  values  for  various  plane  strain,  boundary  value  problems  and  the  extent  of  the  J- 
dominant  region  which  is  synonymous  with  the  HRR  field.  Also  higher  order  terms  for  the 
asymptotic  stress,  strain  and  in  some  cases  the  deformation  fields  for  a  plane  crack  tip 
field  in  a  power  hardening  material  were  derived  by  Li  and  Wang  (1986),  Sharmar  and 
Aravas  (1991),  Yang  et  al  (1993)  and  Zia  et  al  (1993).  Noteworthy  in  these  theoretical 
analysis  is  the  two-parameter,  ductile  fracture  criterion,  which  was  suggested  by  Li  and 
Wang  (1986),  based  on  a  critical  Jic  and  a  high  triaxial  stress  ahead  of  a  plane  strain 
crack  front.  These  two-parameters,  together  with  the  J-integral,  formed  the  basis  of  the 
J-Q  theory  of  fracture  (O'Dowd  and  Shih,  1991,  1992).  More  recently,  Wang  and  Shih 
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(1993)  have  shown  that  the  cumulative  sum  of  the  higher  order  terms  is  equivalent  to  the 
triaxial  stress  and  thus  O'Oowd  and  Shih  (1993)  have  redefined  Q  to  represent  the  entire 
difference  between  the  actual  crack  tip  stress  and  the  HRR  singular  component.  For  small 
scale  yielding,  the  triaxiality  parameter,  Q,  was  related  to  the  elastic  T-stress  and  that 
the  Jic  versus  Q  relation  predicted  from  the  J-Q  ductile  fracture  criterion  correctly 
predicted  the  experimental  results  by  Kirk  et  al  (1993). 

Most  of  the  literature  discussed  above  were  related  to  the  state  of  plane  strain  of  a 
stationary  crack  tip.  The  experimental  investigation,  which  was  described  previously, 
involved  a  near  plane  stress  state  with  an  inevitable  stable  crack  growth.  This  state  of 
plane  stress  effectively  eliminated  the  triaxiality  constraint,  which  is  the  basis  of  the  J- 
Q  theory,  and  thus  the  J-Q  theory  is  not  applicable  for  evaluating  the  test  results  of  thin 
fracture  specimens. 

A  theoretical  analysis  of  the  higher  order  asymptotic  crack  tip  field  for  a  plane 
stress  configuration  has  been  conducted  by  Yang  et  al  (1993a,  1993b).  For  a  power 
hardening  coefficient  of  n  >  3.  their  three  term  displacement  solution  predicted  the  u- 
dispiacement  while  the  v-displacement  required  only  the  HRR  singular  component.  These 
results  are  in  qualitative  agreement  with  the  results  of  the  authors  and  their  colleagues 
with  the  exception  that  the  HRR  displacement  underpredicted  the  measured  u- 
displacement.  In  a  recent  note,  Chao  (1993a)  has  indicated  that  for  a  mode  I  plane  stress 
state  with  n  >  3.2,  the  higher  order  terms  are  uniquely  governed  by  the  J-integral  and  thus 
ductile  fracture  for  such  case  can  be  characterized  by  a  single  parameter  of  J.  This 
important  conclusion  is  in  qualitative  agreement  with  the  results  of  Dadkhah  and 
Kobayashi  (1993)  who  showed  that  throughout  a  stable  crack  growth  of  Aa  ■  2.5  mm,  the 
v-displacements,  which  was  uniquely  represented  by  the  J-integral  values,  in  the  2024-0 
aluminum  alloy,  single  edge  notched  (SEN)  specimens  and  the  biaxially  and  uniaxiaily 
loaded  cruciform  specimens.  Moreover  Chao  (1993)  has  shown  that  the  u-component  of 
the  crack  tip  displacement  approaches  the  HRR  field  with  increasing  biaxiality.  Again, 
this  theoretical  analysis  is  in  qualitative  agreement  with  the  experimental  results  by 
Dadkhah  and  Kobayashi. 

The  above  brief  review  indicates  that  further  theoretical  and  experimental  analyses 
are  needed  before  the  two-parameter  fracture  criterion  such  as  the  J-Q  or  J-T  theory  can 
be  applied  to  a  plane  stress  fracture  specimen.  The  purpose  of  this  paper  is  to  provide 
such  additional  experimental  evidence. 
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TWO-PARAMETER  CRACK  TIP  STATE 

The  high  order  asymptotic  elastic-plastic  crack  tip  stress  field  by  Li  and  Wang 
(1968),  Sharma  and  Aravas  (1991)  and  Yang  et  al  (1993)  is  based  on  the  J2  deformation 
theory  with  the  following  power  hardening  material  of  Ramberg-Osgood; 
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where,  j  -  1  or  2  corresponds  to  a  Cartesian  coordinate  system  with  axes  parallel  or 
perpendicular  to  the  crack,  respectively;  E  and  v  are  the  modulus  of  elasticity  and 
Poisson's  ratio,  respectively;  Oo  and  eo  are  the  yield  stress  and  strain,  respectively;  a  and 
n  are  material  constants;  and  oe.  and  Sjj  are  the  equivalent  and  deviatoric  stress, 
respectively. 

For  a  plane  stress  state.  Gang  et  al  (1994)  has  shown  that  the  two  term 
representation  of  the  asymptotic  crack  tip  field,  in  polar  coordinates,  can  be  represented 
at  the  crack  tip  as 
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i,j  -  1.2  corresponds  to  a  Cartesian  coordinate  system  with  its  origin  at  the  crack  tip;  nj 
represents  a  unit  vector  normal  to  any  contour  which  encloses  the  crack  tip,  u  n,  Uizieiji, 
ei]2 .  oiji  and  oija  are  tabulated  angular  functions  (Yang  et  al,  1993;  Shih,  1983)  of  the 
polar  coordinate,  e,  and  the  hardening  exponent,  n,  and  In.  is  a  tabulated  constant  (Shih 
1983). 

Equations  (1),  (3)  and  (4)  can  be  used  to  evaluate  directly  the  J-integral  from  the 
measured  displacement  field  along  rectangular  contours  which  encompass  the  aack  tip. 
First  the  strain  field  at  each  loading  stage  is  computed  from  the  measured  displacement 
field  at  a  given  point  on  the  contour.  Then  the  corresponding  stress  is  computed  from 
equation  (1)  and  together  with  the  measured  strain,  the  strain  energy  density  is 
determined.  Details  of  this  analysis,  which  is  based  on  the  deformation  theory  of 
plasticity  and  which  was  also  the  basis  of  the  original  definition  of  J-integral  (Rice, 
1968),  and  the  numerical  Integration  procedure  is  given  in  Kang  and  Kobayashi  (1988). 

The  J-integral  value  determined  above  was  then  substituted  into  equation  (2)  to 
compute  the  HRR  singular  displacement  terms.  The  second  order  term  could  not  be 
computed  since  the  corresponding  Ui2  was  not  given  in  Yang  et  al  (1993a).  Assuming  that 
the  second  order  term  is  much  larger  than  the  sum  of  the  remaining  higher  order  terms,  an 
indirect  check  on  the  second  order  term  can  be  made  by  noting  that  it  varies  as  rn/(n>i). 

EXPERIMENTAL  ANALYSIS 

Experimental  Procedure 

The  orthogonal  displacement  components  surrounding  the  crack  tip  were  measured 
via  moir6  interferometry  using  two  specimens,  each  with  a  different  grating.  The  first 
specimen  was  coated  with  a  fine,  cross  diffraction  grating  of  1200  lines/mm  and  was 
used  in  the  initial  phase  of  loading  prior  to  stable  crack  growth.  The  second  specimen 
was  coated  with  a  relatively  coarse,  cross  diffraction  grating  of  40  lines/mm  and  was 
used  to  record  the  crack  displacement  after  stable  crack  growth  commenced.  This  coarse 
grating  was  necessary  due  to  the  gross  yielding  and  the  large  strain  components,  which 
generated  a  moir6  fringe  pattern  too  dense  to  resolve.,  associated  with  the  large  stable 
crack  growth  in  a  ductile  specimen.  The  fine,  cross  diffraction  grating  was  transferred 
using  the  procedure  of  Post  (1993).  The  coarse,  crossed  diffraction  grating  was 
transferred  onto  the  specimen  surface  using  photoresist  and  is  similar  to  the  procedure 
developed  by  Ifju  and  Post  (1991).  However,  in  this  study,  the  highly  polished  surface  of 
the  aluminum  specimen  provided  sufficient  reflectivity  and  thus  an  evaporated  aluminized 
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coating  was  not  used.  This  reflective  specimen  surface  also  eliminated  the  loss  of  moird 
fringes  at  high  strain  where  an  aluminized  coating  may  craze  and  obliterate  the 
diffraction  grating. 

The  specimen  was  illuminated  by  a  four-beam  moir6  interferometer  for 
simultaneous  recording  of  the  two  orthogonal  displacement  fields  (Guo  and  Kobayashi, 
1993).  A  special  four-beam  moir6  interferometry  setup  for  low  frequency  grating  (Wang 
et  al,  1993),  which  can  record  the  moir§  fringes  in  a  severely  warped  specimen  surface, 
such  as  that  at  the  crack  tip,  was  used.  As  shown  in  Figure  1,  moir6  interferometry 
fringes  are  generated  by  the  two  coherent  beams  which  are  projected  on  to  the  deep 
diffraction  grating  on  the  specimen  surface  and  interfere  with  the  radiated  diffracted 
beam.  These  fringes  are  visible  from  all  directions  and  thus  a  good  contrast  of  the  moird 
fringes  in  the  crack  tip  plastic  zone  is  obtained  even  when  the  specimen  is  warped.  Figure 
1  shows  the  moird  interferometry  setup  and  die  u-v  mirror  arrangement  used  in  this 
study.  The  coarse  diffraction  grating  reduced  the  incident  angle  of  the  four  beams  thus 
simplifying  the  u-v  mirror  supports. 

Roeeimen 

The  specimen  consisted  of  a  fatigue  precracked,  thin  single-edged  notch  (SEN), 
2024-T3  aluminum  alloy  specimen  shown  in  Figure  2.  The  moir6  diffraction  grating 
covered  a  region  of  25.4  x  50.8  mm  surrounding  the  crack  as  shown.  Also  shown 
schematically  are  the  integration  contours.  Figure  3  shows  the  uniaxial  stress-strain 
relations  for  two  directions,  i.e.  parallel  and  perpendicular  to  the  rolling  direction  of  the 
sheet.  Also  shown  is  the  average  of  the  two  relations  and  the  two  power  hardening 
coefficients  and  the  least  square  fitted  o-e  relation  (Dadkhah  and  Kobayashi,  1990).  The 
SEN  specimen  was  subjected  to  uniaxial  tensile  loading  in  a  displacement  controlled 
testing  machine  and  the  moird  interferometry  patterns  were  recorded  at  various  stages  of 
stable  crack  growth. 

NUMERICAL  ANALYSIS 

A  commercial  finite  element  code  was  used  to  compute  the  elastic-plastic  state 
associated  with  the  stable  crack  growth  in  this  specimen.  The  objective  of  this 
numerical  analysis  was  to  generate  results,  which  can  be  compared  with  the  experimental 
results,  at  various  stages  of  stable  crack  growth  with  large  scale  yielding. 
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Unlike  the  traditional  finite  element  (FE)  analysis,  the  measured  displacements  near 
the  boundary  of  the  moird  grating  were  used  as  input  boundary  conditions  to  the  FE  model 
of  the  SEN  specimen.  As  discussed  by  Hareesh  and  Chiang  (1989)  and  Sivaneri  et  ai 
(1991),  this  procedure  not  only  resulted  in  saved  computer  time  but  provided  detailed 
information,  which  is  lacking  in  the  moird  analysis  particularly  at  the  initial  stage  of 
stable  crack  growth,  in  the  immediate  vicinity  of  the  crack  tip.  Also,  as  described  in 
Hareesh  and  Chiang  (1989),  this  ioad  path  dependent,  elastic-plastic  finite  element 
analysis  must  start  from  an  early  stage  of  plastic  yielding  and  proceed  incrementally 
along  the  loading  path. 

Figure  4  shows  the  finite  element  model  used  in  this  analysis.  The  power  hardening 
stress-strain  relation  used  in  this  incremental  elastic-plastic  analysis  was  that  of 
Dadkhah.  Due  to  the  sensitivity  of  the  FE  method  to  displacement-prescribed  boundary 
conditions,  a  second  order  curve,  which  was  fitted  to  the  measured  boundary 
displacements  obtained  from  moird  interferometry,  was  used  as  the  input  boundary 
condition. 

RESULTS 

Eleven  increments  of  load,  i.e.  five  increments  for  the  specimen  with  a  fine  moir^ 
gratings  and  six  increments  for  the  specimen  with  a  coarse  grating,  were  applied  and  the 
corresponding  moird  interferometry  fringes  were  recorded.  Figure  5  shows  the  measured 
load  versus  the  computed  load  line  displacement  for  the  two  specimens  with  the  fine  and 
the  coarse  grating.  The  continuity  in  the  two  relations  Justify  the  use  of  the  two 
specimens  in  this  study.  The  J  values  were  then  computed  along  three  contours  of  10mm  x 
10mm,  15mm  x  15mm  and  17.5mm  x  25mm  for  each  increment  of  loading. 

Figures  6a  and  6b  show  typical  moird  fringe  patterns  corresponding  to  the 
displacement  parallel,  u,  and  perpendicular,  v,  to  the  crack  using  the  fine  and  coarse 
diffraction  gratings,  respectively.  Also  shown  in  Figs.  6a  and  6b  are  the  rectangular 
contours  used  for  the  J-integral  computation.  This  computation  for  the  specimen  with  a 
coarse  diffraction  grating,  using  moir8  interferometry  data,  was  conducted  during  the 
later  stage  of  loading  which  yielded  denser  moir4  fringe  patterns. 

Figures  7a,  7b  find  7c  show  a  comparison  of  the  various  J  values  of  interest.  Jfem 
designates  the  value  of  the  J-integral  computed  by  the  commercial  FE  code,  using  the 
moir6  fringe  data  as  input  boundary  conditions,  prior  to  and  with  crack  extension.  Jlefm 
was  computed  by  J-G  (strain  energy  release  rate)  based  on  the  K|  from  the  stress 
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intensity  factor  listing  in  Tada  et  al  (1973).  Jshih  was  obtained  using  the  procedure  and 
tables  of  Kumar  et  al  (1981).  J  values  were  also  obtained  using  an  algorithm  which 
calculated  the  value  of  the  J-integral  along  a  designated  integration  contour  on  the  moir^ 
fringe  pattern.  Jexperimental  is  the  value  of  J  associated  with  the  largest  contour.  Figure 
7a  shows  that  Jexperimental  is  in  good  agreement  with  Jlefm  at  lower  loading  as  expected 
and  increasingly  differs  at  a  higher  loading  prior  to  stable  crack  growth.  Figure  7b  shows 
that  the  value  of  Jexperimental  is  close  to  the  value  of  Jshih  at  higher  loads  and  during 
stable  aack  growth.  Another  point  of  interest  is  that  Jexperimental  becomes  highly 
contour  dependent  before  a  crack  extension  of  aa  >  1.5mm  is  reached  as  shown  in  Fig.  7c. 
This  is  unexpected  in  light  of  previous  numerical  studies  by  Shih  et  al  (1981). 

Figures  8a  and  8b  show  the  log-iog  plots  of  crack  tip  u-  and  v-dispiacements, 
respectively,  prior  to  stable  crack  growth  for  three  angular  orientations.  While  the 
measured  and  LEFM  u-  and  v-displacements  are  in  reasonable  agreement,  as  expected  at 
this  low  load,  a  similar  agreement  is  not  observed  in  the  very  vicinity  of  the  crack  tip. 

Due  to  the  uncertainty  in  the  rigid  body  displacement  component  in  the  measured  u- 
displacement,  only  the  slope  and  not  the  absolute  value  of  this  displacement  can  be 
compared  with  the  slope  of  the  LEFM  component,  i.e.  0.5,  and  with  the  slope  of  the  HRR 
component  which  is  0.083. 

Figures  9a  and  9b  show  the  variations  in  the  two  normal  strains,  exx  and  eyy,  with 
crack  tip  distance,  r,  at  three  angular  orientations  of  e  -  30°,  45°  and  60°  prior  to  stable 
crack  growth.  While  the  eyy  followed  the  LEFM  eyy  as  expected  from  the  results  of  Fig.  8b, 
the  exx  deviated  from  the  LEFM  exx  at  r  <  1  mm  in  the  region  of  e  <  45°-  Such  difference 
was  unexpected  since  Fig.  8a  showed  that  the  slope  of  the  log-log  plot  of  the  u- 
displacement  was  the  anticipated  LEFM  slope  of  1/2. 

Figures  10a  and  10b  show  the  log-log  plots  of  the  crack  tip  u-  and  v-displacements, 
respectively,  at  the  onset  of  stable  crack  growth.  The  u-displacement  continues  to 
maintain  its  elastic  response  while  the  v-displacement  is  approaching  the  HRR 
displacement  at  about  r  -  1  mm. 

Figures  I1a  and  11b  show  the  variations  in  the  two  normal  strains,  exx  and  cyy,  with 
r  at  e  -  30°,  45°  and  60°.  Also  plotted,  in  addition  to  the  corresponding  LEFM  normal 
strain,  are  the  HRR  strains  computed  using  the  measured  J-integral  values.  Despite  the 
increasing  size  in  the  plastic  region,  these  results  are  similar  to  those  in  Figs  9a  and  9b 
where  the  LEFM  strain  field  still  dominates. 


Figures  12a  and  12b,  as  well  as  13a  and  13b,  show  the  log-log  plots  of  the  crack  tip 
displacements  for  a  moderate  stable  crack  growth,  i.e.  Aa  «  1.5  mm  and  3.9  mm, 
respectively.  The  slope  and  the  value  of  the  v-displacement  is  now  in  good  agreement 
with  those  of  the  HRR  displacement.  The  slope  and  the  value  of  u-displacement,  however, 
continues  to  follow  those  and  the  value  of  the  LEFM  displacement. 

The  Q  value  in  these  figures  represents  the  absolute  value  of  the  difference  between 
the  measured  and  the  HRR  displacements  computed  from  the  measured  J  values  and  is  not 
the  Q  value  used  by  Shih  et  al.  The  simplified  J-Q  theory  predicts  for  a  state  of  plane 
stress,  vanishing  Q  components  as  per  Aravas  (1993).  Figures  12  and  13  show  Q 
components  of  the  same  order  as  the  measured  u-  and  v-  values  The  slopes  of  the  Q 
components  for  u-  and  v-displacements  vary  from  0  to  0.18  within  the  region  of  roo/J  <  3. 
In  the  region  of  3  <  roo/J  <10,  the  average  slope  of  the  O  components  for  the  u-  and  v- 
displacements  are  approximately  0.9  and  0,  respectively.  While  these  values  are  not 
exactly  the  n/(n-i-l)  -  0.92  predicted  by  Gang  et  al  (1994),  a  qualitative  agreement  is 
noted. 

Figures  14a  and  14b  and  Figs.  15a  and  15b  show  the  variations  in  two  crack  normal 
strains,  exx  and  cyy,  for  moderate  crack  growths  of  Aa  -  1.5  and  3.9  mm,  respectively.  Also 
plotted  are  the  corresponding  HRR  strains  whit^  were  computed  by  using  the 
experimentally  determined  J-integral  values.  In  terms  of  absolute  values,  the  HRR  strain 
grossly  over  predicts  exx  while  it  agrees  well  with  the  experimental  eyy  for  e  ■>  60°  but 
not  so  well  at  e  -  30°  and  45°. 

Figure  16  shows  the  variations  in  the  crack  tip  opening  angle  (CTOA),  as  measured  1 
mm  from  the  crack  tip,  using  the  coarse  and  fine  grating  moire  interferometry  data,  with 
the  increase  in  loadline  displacement.  The  CTOA  reached  a  nearly  constant  value  at  a 
crack  extension  of  Aa  «  1.5  mm  with  a  leveled  resistance  curve  thereafter.  Since  CTOA  is 
related  to  the  maximum  strain  at  the  crack  tip,  this  result  suggests  that  CTOA  could  be  a 
viable  ductile  fracture  parameter.  The  well-known  disadvantage  is  that  CTOA  cannot  be 
readily  measured  with  a  normal  testing  facility  and  that  it  is  impractical  to  determine  in 
complex  flaw  geometries,  such  as  a  surface  flaw. 

CONCLUSIONS 

Limited  experimental  results  involving  the  crack  tip  displacement  fields  in  thin 
aluminum  SEN  with  small  stable  crack  growth  showed  that  the  J-Q  theory,  based  on  the 
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plane  strain  form  of  O'Dowd  and  Shih  (1993),  may  not  be  present.  On  the  other  hand,  the 
crack  tip  field  of  Yang  et  al  (1993)  could  exist. 

For  an  elastic-plastic  fracture  analysis  of  thin  plates  in  the  presence  of  stable 
crack  growth,  the  present  J-integral  computation  procedures  based  on  Kumar  et  al  (1981) 
and  on  a  commercial  code  must  be  reinvestigated  in  view  of  the  observed  large 
discrepancies. 

Preliminary  results  show  that  the  CTOA  could  be  a  viable  ductile  fracture  parameter 
for  stable  crack  growth  and  ductile  fracture. 
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Figure  2.  2024-T3  aluminum  alloy  SEN  specimen  with  J-integral  contours  shown 
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Figure  3.  Uniaxial  stress-strain  curve  for  2024-T3 


Rgure  4.  Finite  element  mesh  used  in  elastic-plastic  analysis 
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Rgure  6a.  Moire  interferometry  patterns  of  the  2024-T3  SEN  specimen  with 
J  contours  shown  [Load«689.5  (N)].  Fine  grating  analysis 
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Figure  6b.  Moire  interferometry  patterns  of  the  2024-T3  SEN  specimen  with 
J  contours  shown  [Loadsf  134  (N)].  Coarse  grating  analysis 


Variation  in  the  J-integral  values  of  2024-T3  specimen  under  stable  crack  growth 


Figure  8a.  Crack  tip  u-dlsplacement  of  2024-T3  SEN  specimen  prior  to  stable  crack  growth.  Fine  grating  analysis 


Rgure  8b.  Crack  tip  v-displacement  of  2024-T3  SEN  spedmen  prior  to  stable  crack  growth.  Firre  grating  analysis 


EXX 


roq/J 

0  200  400  600  800  1000  1200 


rao/J 

xlO-9  0  100  200  300  400  SOO 


r  (mm) 


xlO-3 


rcJo^J 

50  too 


150 


-1  h 


■2  h 


£xx 


□  □ 


EXPERIMENTAL  STRAINS^ 
LEFM  STRAINS 
LOADS  690  (N) 

Aas  0  (mm) 

J0xp^35  (PaHn) 

es60o 


-5  h 


I  -  I  .  I  .  I  ■  I  ■  I 
02  0.4  02  0.8  1.0  12  1.4 

r(mm) 


Figure  9a.  Crack  tip  Exx  of  2024-T3  SEN  specimen  prior  to  stable 
crack  growth.  Fine  grating  analysis 
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Figure  db.  Crack  tip  £yy  of  2024-T3  SEN  specimen  prior  to 
stable  crack  growth.  Rne  grating  analysis 
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Rgure  1 1  a.  Crack  tip  Exx  of  2024-T3  SEN  specimen  at  the  onset  of 
stable  crack  growth.  Rne  grating  analysis 
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Figure  1 1  b.  Crack  tip  £yy  of  2024-T3  SEN  specimen  at  the  onset  of 
stable  crack  growth.  Fine  grating  analysis 


Figure  12b.  Crack  tip  v-displacement  of  2024-T3  SEN  specimen  with  smati  crack  growth.  Coarse  grating  analysis 
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Rgure  13b.  Crack  tip  v-dlsplacement  of  2024-T3  SEN  specimen  with  moderate  stable  crack  growth. 
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Figure  1 4a.  Crack  tip  Exx  of  2024-T3  SEN  specimen  with  small 
stable  crack  growth.  Coarse  grating  analysis 


Rgure  14b.  Crack  tip  eyy  of  2024-T3  SEN  specimen  with  small 
stable  crack  growth.  Coarse  grating  analysis 
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Figure  15a.  Crack  tip  Exx  of  2024-13  SEN  specimen  with  moderate 
stable  crack  growth.  Coarse  grating  analysis 
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Rgure  15b.  Crack  tip  Eyy  of  2024-T3  SEN  specimen  with  moderate 
stable  crack  growth.  Coarse  grating  analysis 
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